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ABSTRACT 

The origin of the extragalactic gamma-ray background radiation at 1-10 MeV is still unknown. 
Although the cosmic X-ray background (CXB) up to a few hundreds keV can be accounted for by the 
sum of Active Galactic Nuclei (AGNs) , current models of AGN spectra cannot explain the background 
spectrum beyond ~1 MeV, because of the thermal exponential cutoff of electron energy distribution 
assumed in the models. Here we construct a new spectral model by calculating the Comptonization 
process including nonthermal electrons, which are expected to exist in an AGN hot corona if it is heated 
by magnetic reconnections. We show that the MeV background spectrum can nicely be explained by 
our model, when coronal electrons have a nonthermal power-law component whose total energy is a 
few percent of the thermal component and whose spectral index is din N e /d\nE e w —4. Although 
the MeV gamma-ray flux from such a component in nearby AGN spectra is below the detection limit 
of past observations, it could be detected by planned future MeV detectors. We point out that the 
amount of the nonthermal component and its electron index are similar to those found for electrons 
accelerated by magnetic reconnections in solar flares and the Earth magnetosphere, giving a support 
to the reconnection hypothesis for the origin of hot AGN coronae. 
Subject headings: diffuse radiation — galaxies: active — gamma rays: theory 



1. INTRODUCTION 

It is well-known that normal active galactic nuclei 
(AGNs) explain the cosmic X-ray background (CXB) be- 
low several hundreds keV (for reviews see Boldt 1987; 
Fabian & Barcons 1992; Ueda et al. 2003, hereafter U03; 
Gilli, Comastri, & Hasinger 2007). It is also known that 
rare AGNs of the blazar type make a considerable con- 
tribution to the cosmic gamma-ray background in the 
energy range from 100 MeV to 100 GeV, which has a 
hard power-law spectrum (almost flat in the vF v plot), 
though blazars may not explain all of the background 
flux, leaving some room for possible contributions from 
completely different sources (Narumoto & Totani 2006, 
and references therein). 

The origin of the gamma-ray background at the gap be- 
tween these two energy regions, i.e., ~1-10 MeV, has also 
been an intriguing mystery. The AGN spectra adopted in 
population synthesis models of the CXB cannot explain 
this component because of the assumed exponential cut- 
off at a few hundred keV. The background spectrum in 
the 1-10 MeV band is much softer (photon index a ~2.8) 
than the GeV component, indicating a different origin 
from that above 100 MeV (e.g. Sreekumar et al. 1998). 

A few candidates have been proposed to explain the 1- 
10 MeV background. One is the nuclear decay gamma- 
rays from type la supernovae (SNe la) (Clayton & Ward 
1975; Zdziarski 1996; Watanabe et al. 1999). However, 
recent studies based on the latest measurements of the 
cosmic SN la rates show that the background flux ex- 
pected from SNe la is about an order of magnitude lower 
than observed (Ahn, Komatsu, & Hoflich 2005; Strigari 
et al. 2005). There is a class of blazars called "MeV 
blazars", whose spectra have peaks at ~ MeV (Blom 
et al. 1995; Sambruna et al. 2006), and these MeV 
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blazars could potentially contribute to the MeV back- 
ground. However, quantitative estimate of the contri- 
bution is difficult because of the still small sample avail- 
able at present. Annihilation of the dark matter particles 
has also been discussed (Ahn & Komatsu 2005a, 2005b; 
Rasera et al. 2006; Lawson & Zhitnitsky 2007), but there 
is no natural particle physics candidate for such a dark 
matter particle with a mass scale of ~ MeV. The mo- 
tivation of MeV dark matter has been inspired by the 
511 keV line emission from the Galactic center, but a 
few astrophysical explanations are possible for this line 
emission (Totani 2006, and references therein). 

An important feature of the MeV background spec- 
trum is that its power-law spectrum is smoothly con- 
nected to the peak of the CXB spectrum. If the origin 
of the MeV background is completely different from that 
of the CXB, such a smooth connection would be sur- 
prising. Rather, it seems more plausible that the MeV 
background flux is composed of the same populations 
that make the CXB, and simply the current AGN spec- 
tral models are not sufficient to describe the MeV spec- 
tra. The X-ray AGN spectra are well described by the 
Comptonization of seed UV photons by the hot coronal 
electrons (Zdziarski et al. 1994, 1995), and the cut-off 
at ^100 keV reflects the thermal energy distribution of 
the hot electrons. Although the AGN spectra indeed 
show evidence for such a cut-off (Zdziarski et al. 1995; 
Zdziarski, Poutanen, & Johnson 2000), a small amount 
of additional non-thermal electrons with a soft spectrum 
is sufficient to explain the MeV background. Due to the 
limited sensitivity of current MeV gamma-ray observa- 
tions, the presence of such non-thermal components is 
not strongly constrained even in the spectra of nearby 
brightest AGNs. Furthermore, it is believed that coro- 
nae around accretion disks share some common features 
with the solar corona (e.g., Galeev, Rosner, & Vaiana 
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1979), and magnetic reconnection in AGN coronae is a 
good candidate for the origin of hot electrons (Liu, Mi- 
neshige, & Shibata 2002). It is well known that particles 
are accelerated to nonthermal energies by reconnections 
in solar flares (Shibata et al. 1995). 

Here we construct a new model of the X/gamma-ray 
spectra of AGNs, by calculating the Comptonization pro- 
cess by hot electrons having both thermal and nonther- 
mal components. We also calculate the CXB spectrum 
based on our model with the latest knowledge of the 
cosmological evolution of the AGN luminosity function, 
and determine the amount and spectrum of the nonther- 
mal electrons in AGN coronae to explain the MeV back- 
ground. We discuss the implied nature of nonthermal 
electrons in the context of the reconnection heating sce- 
nario of the AGN coronae, comparing our results with 
those found in the reconnections occurring in the solar 
flares and the Earth magnetosphere. 

Rogers & Field (1991) and Field & Rogers (1993) pre- 
sented an AGN spectral model that can explain the 
MeV background spectrum by nonthermal relativistic 
electrons. However, their model only considers the non- 
thermal component without a thermal component, and 
it requires a lower cut-off of j e ~ 30 in the nonther- 
mal component, which is difficult to interpret as they 
mentioned in their paper. Our model considers both the 
thermal and nonthermal coronal electrons whose spectra 
are smoothly connected to each other, which is a natu- 
ral extension of the popular AGN spectral models in the 
recent literature. Stecker, Salamon, & Done (1999) also 
discussed a possibility that the MeV background is ex- 
plained by nonthermal tails in AGN spectra, quoting the 
spectrum of the Galactic stellar-mass black hole candi- 
date Cyg X-l. However, a physical model to explain the 
nonthermal tail in an AGN spectrum was not presented. 

Throughout this paper, we adopt the cosmological pa- 
rameters of (h , ft m , ^a) = (0. 7,0.3,0. 7). 

2. MODEL DESCRIPTION 

2.1. the AGN Spectra with Nonthermal Coronal 
Electrons 

The main shape of X-ray AGN spectra is deter- 
mined by Comptonization of UV photons emitted from 
optically-thick accretion disks by hot electrons in coro- 
nae. As in many previous studies, we consider a simple 
spherical and uniform distribution of the coronal elec- 
trons. The seed UV photons are injected at the cen- 
ter and then become X-ray photons when they escape 
the coronal region after Comptonization. In addition to 
the hot thermal electrons assumed in the conventional 
X-ray spectral models of AGNs, we introduce higher en- 
ergy nonthermal electrons in AGN coronae, whose en- 
ergy distribution is a power-law as dN e /dE e oc E~ T '. We 
introduce the transition electron Lorentz factor 7 tr , cor- 
responding to the transition electron energy E e = rn e 7t r 
where the electron spectrum dN e /dE e has the same value 
for the thermal and nonthermal components. This j tT is 
the lower limit of the Lorentz factor distribution of the 
nonthermal component, and hence there are no nonther- 
mal electrons at E e < m e 7 tr . We also set an upper bound 
as j u = 10 5 , although this hardly affects our results if the 
maximum photon energy well extends beyond 10 MeV. 

We set the coronal temperature to be kT e = 256 keV 
and assume a blackbody spectrum for UV seed photons 



from a cooler disk with Td = 10 eV, following the con- 
ventional thermal models (e.g. Zdziarski et al. 1994). 
The degree of Comptonization is determined by the op- 
tical depth for Thomson scattering, tt- We found that 
the spectral photon index in the X-ray band, ax, be- 
comes close to that typically found in observed spectra 
(ax ~ L9, e.g., Nandra & Pounds 1994; Turner et al. 
1997; George et al. 1998) when we set r T = 0.24. There- 
fore we use this value throughout this letter; this value 
is also similar to those used in the conventional mod- 
els. It should be noted that ax is hardly changed even 
if we introduce the nonthermal electron component with 
an amount that is necessary to explain the cosmic MeV 
background. 

We then trace the Comptonization process using a 
Monte Carlo method. The calculation method used here 
is mainly based on that in Pozdniakov, Sobol, & Siuni- 
aev (1977) and Gorecki & Wilczewski (1984), but their 
original formalism in the laboratory frame is not opti- 
mized for the ultra-relativistic region. To calculate the 
scattering by high energy nonthermal electrons more ef- 
ficiently, we added a new formulation in the rest frame 
of relativistic electrons based on Corman (1970). 

The reflection of X-ray photons by cool, optically thick 
matter is also an important feature of AGN X-ray spec- 
tra. We calculate this by using the PEXRAV model 
(Magdziarz & Zdziarski 1995) in the XSPEC package as 
done in U03. Because of the limitation for the acceptable 
input spectrum in PEXRAV, we use a power-law spec- 
trum (ax = 1.9) plus an exponential cutoff at E = 500 
keV, which is a good approximation of the Comptonized 
spectrum only with the thermal electrons. The newly 
added nonthermal electrons would change the spectrum 
significantly only at E > 1 MeV, and hence this treat- 
ment is appropriate for the reflection component which 
is important only at ~ 1-100 keV. 

2.2. Calculation of the Cosmic Background Radiation 

We calculate the CXB spectrum by integrating our 
AGN spectral model in the redshift and luminosity space, 
using the X-ray AGN luminosity function of U03. Fol- 
lowing the same formulation given in U03, we take 
into account the absorption column density distribution 
(Nh function) and the contribution from Compton-thick 
AGNs. We confirm that our main conclusion hardly 
change if instead we use a more recent population syn- 
thesis model by Gilli et al. (2007), as described below. 

3. RESULTS 

Figure [T] shows the models of AGN spectra calcu- 
lated according to the procedures in the previous sec- 
tion. Here, we do not take into account the reflec- 
tion component and the absorption effect, to show the 
pure spectrum of the Comptonization. We set T = 3.8 
and 7t r = 4.4 as our standard model (solid line), be- 
cause we will find that these values give the best-fit 
MeV background spectrum to the data. In this standard 
model, 3.5% of the total electron energy is carried by the 
nonthermal electrons. To illustrate effects of changing 
parameters, we also show the spectra with parameters 
slightly changed from those for the best fit model. The 
conventional case only with the thermal electrons is also 
shown, where an exponential cutoff above ~100 keV is 
seen as expected. 
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Fig. 1. — The AGN spectra in X-ray and gamma-ray bands 
calculated by our model. The flux is shown in an arbitrary unit of 
E 2 dN/dE, where dN/dE is a differential photon spectrum. They 
are Comptonization of UV seed photons without taking into ac- 
count the reflection component and the absorption effect. The 
thick solid curve is our standard spectrum with T = 3.8 and 
7t r = 4.4. The other thick curves are for the cases of different model 
parameters as indicated in the figure. The thick dotted curve is the 
spectrum only with the thermal component (kT e = 256 keV). The 
thin dotted curve is the input UV spectrum (a black body with 
T d = 10 eV). 

Figure [2] shows the cosmic background radiation in- 
tegrated with the luminosity function of U03 using our 
AGN spectral model. It has been known that the pre- 
dicted CXB spectrum below 100 keV by the U03 model 
is 10-20% higher than the HEAO-1 data. The origins 
of this discrepancy are still controversial. The intensity 
and shape of the CXB in the 10 keV - 1 MeV band could 
be uncertain at «20% level 1 . The population synthe- 
sis models have also uncertainties, such as the intensity 
of the reflection component assumed in the AGN spec- 
tra, the (unknown) number density of Compton-thick 
AGNs, and the parameters of the luminosity function 
itself that could be subject to cosmic variance in deep 
surveys. However, here we emphasize that the uncer- 
tainty hardly affects our conclusions. To confirm this, 
we also calculate the MeV gamma-ray background flux 
with the U03 luminosity function by artificially lowering 
its normalization by 20% to match the HEAO-1 data. 
We find that the best fit value of 7t r is slightly changed 
by 0.4, while the change of T is negligible. 

As mentioned above, we find here that the cosmic back- 
ground spectrum from X-ray to 10 MeV band can nicely 
be explained by our model with T = 3.8 and 7 tr = 4.4. 
The next question is then how natural are these param- 
eters in the context of the theoretical picture of hot elec- 
trons in AGN coronae. We will discuss about this issue 
in the next section. 

4. DISCUSSIONS 

4.1. Reconnection and Nonthermal Electrons 

As discussed in §1, magnetic reconnection is the pri- 
mary candidate for the origin of the nonthermal elec- 
trons in AGN coronae, and we can compare the inferred 
amount and spectrum of electrons with those observed in 

1 In Figure 2, a discontinuity is seen between the HEAO-1 A4 
MED data above 100 keV and the HEAO-1 A4 LED data below 
100 keV (Gruber et al. 1999, see also ; Revnivtsev et al. 2005; 
Churazov et al. 2007; Frontera et al. 2007). If we use the model 
by Gilli et al. (2007), the discrepancy is reduced below 100 keV 
but is enhanced above 100 keV. 
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Fig. 2. — The spectrum of the cosmic background radiation 
in X-ray and gamma-ray bands, predicted by our model of AGN 
spectra shown in Fig. [l] For each line-marking, the corresponding 
AGN spectrum in Fig. [T] is used for the calculation. The data 
points of HEAO-1 (Gruber et al. 1999) SMM (Watanabe et al. 
1999), and COMPTEL (Kappadath et al. 1996) experiments are 
also shown. For the SMM data the thin dotted line indicates the 
±ler uncertainty region. 

reconnections in other objects. In solar flares, the elec- 
tron spectrum accelerated by reconnection and injected 
into the solar surface (foot points) can be estimated by 
thick-target bremsstrahlung (Brown 1971; Piana 1994), 
and a value of T ~ 4 is inferred for giant solar flares 
(Holman et al. 2003; Lin 2006). In the reconnection 
diffusion region of the Earth's magnetotail, T = 3.8 has 
been measured by the Wind spacecraft (0ieroset et al. 
2002). Interestingly, these values are very similar to what 
we found to explain the MeV background spectrum by 
AGNs. 

The relative amounts of the thermal and nonthermal 
electrons are difficult to predict, but it should be deter- 
mined by the balance between the cooling of the ther- 
mal electrons, thermalization of nonthermal electrons, 
and energy input rate by reconnections. In solar flares, 
the total energy of accumulated nonthermal electrons is 
comparable with or even larger than that of the thermal 
electrons (Holman et al. 2003). In the directly observed 
electron spectrum in the Earth's magnetotail, the non- 
thermal component is smoothly connected to the ther- 
mal component at the energy where the thermal elec- 
tron spectrum declines by the exponential cutoff, which 
is reminiscent of the cosmic background spectrum from 
X-ray to the MeV band. Therefore, the inferred amount 
of nonthermal electrons in AGN coronae seems quite rea- 
sonable in comparison with the reconnections observed 
in the Sun or the Earth magnetosphere. 

As in the solar flare, some of the nonthermal electrons 
in AGNs would be injected into accretion discs along the 
magnetic field lines, and would produce bremsstrahlung 
emission. Such emission may also contribute to the 
AGN spectra. However, we consider here the inverse- 
Compton scattering since X-ray AGN spectra are well ex- 
plained by the Comptonization process, indicating that 
the efficiency of Comptonization is much higher than 
bremsstrahlung in AGNs, in contrast to the case of X-ray 
radiation from solar flares. It should be noted that the 
radiation efficiency of the thick-target bremsstrahlung 
is generally much lower than the unity because the 
electrons injected into dense material would lose their 
energy by ionization or Coulomb losses rather than 
bremsstrahlung . 



4 



4.2. Implications for Future MeV Observations of AGNs 

Our results suggest that the majority of AGN popula- 
tions that are responsible for the CXB should commonly 
have a nonthermal component beyond the thermal expo- 
nential cutoff in their gamma-ray spectra. Future sensi- 
tive MeV observations may reveal the nonthermal com- 
ponents from nearby AGN spectra. For instance, the 
Advanced Compton Telescope (ACT, see Boggs (2006)), 
scheduled for launch around 2015, will have a continuum 
sensitivity in the 0.2-10 MeV band with E*dF 1 /dE- f ~ 
1 x 10~ 5 (£ 7 /MeV) MeV cm^s" 1 , where dF 1 /dE 1 is the 
differential photon flux. Using the hard X-ray flux of 
NGC 4151 (Sazonov et al. 2007), the brightest Seyfert 
galaxy, we estimate the expected nonthermal MeV flux 
as - 3 x 10- 5 (£ 7 /MeV)-°- 8 MeV cm^s" 1 , which can 
well be detected with the ACT. The detection of such 
components from nearby bright AGNs would give a clear 
test for our model. 

5. CONCLUSIONS 

In this paper, we have shown that MeV gamma-ray 
background can be explained by the same population of 
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